Relativistic static thin dust disks with an inner edge: 
An infinite family of new exact solutions 
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An infinite family of new exact solutions of the Einstein vacuum equations for static and axially 
symmetric spacetimes is presented. All the metric functions of the solutions are explicitly computed 
and the obtained expressions are simply written in terms of oblate spheroidal coordinates. Further- 
more, the solutions are asymptotically flat and regular everywhere, as it is shown by computing all 
the curvature scalars. These solutions describe an infinite family of thin dust disks with a central 
inner edge, whose energy densities are everywhere positive and well behaved, in such a way that 
their energy-momentum tensor are in fully agreement with all the energy conditions. Now, although 
the disks are of infinite extension, all of them have finite mass. The superposition of the first mem- 
ber of this family with a Schwarzschild black hole was presented previously [G. A. Gonzalez and 
A. C. Gutierrez-Pineres, arXiv: 0811.3002vl (2008)], whereas that in a subsequent paper a detailed 
analysis of the corresponding superposition for the full family will be presented. 

PACS numbers: 04.20.-q, 04.20.Jb, 04.40.-b 



I. INTRODUCTION 

The observational data supporting the existence of 
black holes at the nucleus of some galaxies is today so 
abundant, with the strongest dynamical evidence com- 
ing from the center of the Milky Way, that there is no 
doubt about the relevance of the study of binary systems 
composed by a thin disk surrounding a central black hole 
(see 0, for recent reviews on the observational evi- 
dence). Accordingly, a lot of work has been developed 
in the last years in order to obtain a better understand- 
ing of the different aspects involved in the dynamics of 
these systems. Now, due to the presence of a black hole, 
the gravitational fields involved are so strong that the 
proper theoretical framework to analytically study this 
configurations is provided by the general theory of rela- 
tivity. Therefore, a strong effort has been dedicated to 
the obtention of exact solutions of Einstein equations cor- 
responding to thin disklike sources with a central black 
hole (see [3,|j| for thoroughly surveys on the subject). 

Stationary and axially symmetric solutions of the Ein- 
stein equations are of obvious astrophysical importance, 
as they describe the exterior of equilibrium configurations 
of bodies in rotation. At the same time, such spacetimes 
are the best choice to attempt to describe the gravita- 
tional fields of disks around black holes in an exact ana- 
lytical manner. So, through the years, several examples 
of solutions corresponding to black holes or to thin disk- 
like sources has been obtained by many different tech- 
niques. However, due to the nonlinear character of the 
Einstein equations, solutions corresponding to the super- 
position of black holes and thin disks are not so easy to 
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obtain and so, until now, exact stationary solutions have 
not been obtained. 

On the other hand, if we only consider static config- 
urations, the line element is characterized only by two 
metric functions. So, in the vacuum case, the Einstein 
equations implies that one of the metric functions satis- 
fies the Laplace equation whereas that the other one can 
be obtained by quadratures. Furthermore, as the sources 
are infinitesimally thin disks, the matter only enters in 
the form of boundary conditions for the vacuum equa- 
tions. Therefore, as a consequence of the linearity of the 
Laplace equation, solutions corresponding to the super- 
position of thin disks and black holes can be, in principle, 
easily obtained. 

However, if we consider thin disks that extend up to 
the event horizon, the matter located near the black hole 
will moves with supcrluminal velocities, as was shown 
by Lemos and Letelier 0, 0, 0]. So, in order to pre- 
vent the appearance of tachyonic matter, the thin disks 
must have an inner edge with a radius larger than the 
photonic radius of the black hole. Then, the boundary 
value problem for the Laplace equation is mathematically 
more complicated and thus only very few exact solutions 
had been obtained. These kind of solutions were first 
studied by Lemos and Letelier Q by making a Kelvin 
transformation in order to invert the Morgan and Mor- 
gan [8| family of finite thin disks. Now, although the 
second metric function of this solution can not be ana- 
lytically obtained, their main properties were extensively 
analyzed in a series of papers by Semerak, Zacek and 
Zellerin [1 El El El El EI El, by using numerical 
computation when was needed. 

Besides the Lemos and Letelier inverted disks, only 
two other solutions for static thin disks with an in- 
ner edge have been obtained, a first one with inverted 
isochrone disks [l6| and a second one for disks with 
a power-law density [13]. Also, a stationary superpo- 
sition was obtained by Zellerin and Semerak [18| by 
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using the Belinskii-Zakahrov inverse-scattering method, 
but the analysis of its properties is complicated by the 
fact that the metric functions can not be analytically 
computed. Furthermore, this solution involves an un- 
physical supporting surface between the black-hole hori- 
zon and the disk [l9(. Finally, a general class of sta- 
tionary solutions was presented by Klein f20j, by using 
the Riemann-surface techniques, with which physically 
acceptable black hole disk systems can in principle be 
found. 

Now, a common feature of all the above mentioned so- 
lutions is that their metric functions can not be fully ana- 
lytically computed. Thus then, the analysis of their phys- 
ical and mathematical properties is very complicated. 
Moreover, almost all of these solutions do present sin- 
gularities at the inner edge of the disk, perhaps the only 
exemption being the class of solutions presented by Klein 
pol ] . Now then, the ubiquitous presence of this singular- 
ity at almost all the obtained solutions has been consid- 
ered by some authors as tightly connected with the un- 
physical infinite thinness of the source. However, as we 
will show at this paper, if the corresponding boundary 
value problem is properly solved, it is possible to obtain 
a whole family of singularity free solutions. 

In this paper we present an infinite family of new ex- 
act solutions for static thin dust disks with a central in- 
ner edge. These solutions describe disks whose energy 
densities are everywhere positive and well behaved, in 
such a way that their energy-momentum tensor are in 
fully agreement with all the energy conditions. Now, al- 
though the disks are of infinite extension, all of them 
have finite mass. Furthermore, the solutions are asymp- 
totically flat and regular everywhere, as it is shown by 
computing all the curvature scalars. The superposition 
of the first member of this family with a Schwarzschild 
black hole was presented previously [21( , whereas that in 
a subsequent paper a detailed analysis of the correspond- 
ing superposition for the full family will be presented. 

The paper is organized as follows. First, in Sec. 
HI! we present the formulation of the Einstein equa- 
tions for static axially symmetric spacetimes with an in- 
finitesimally thin disk as source. We also present the 
proper boundary conditions and their relationship with 
the physical quantities characterizing the sources. Then, 
in Sec. IIII1 we introduce oblate spheroidal coordinates 
with the ranges choosen in such a way that they be nat- 
urally adapted to the geometry of a thin disk with a cen- 
tral inner edge. The Einstein equations are then solved 
and the metric functions of whole the family of solutions 
are explicitly computed and the obtained expressions are 
simply written in terms of the oblate spheroidal coordi- 
nates. The behavior of the solutions is then analysed in 
Sec. IIVI The analysis is made by studying the behavior 
of all the curvature invariants for the full family as well 
as the behavior of the corresponding energy densities and 
azimuthal pressures. Also we analyze the mass densities 
of the disks and the finitcness of their corresponding total 
mass. Finally, in Sec. El we summarize the results. 



II. THE EINSTEIN EQUATIONS WITH THIN 
DISKLIKE SOURCES 

In order to formulate the Einstein equations for static 
axially symmetric spacetimes with an infinitesimally thin 
disk as source, first we introduce coordinates x a = 
(t, tp, r, z) in which the metric tensor only depends on 
r and z. We assume that these coordinates are quasi- 
cylindrical in the sense that the coordinate r vanishes on 
the axis of symmetry and, for fixed z, increases mono- 
tonically to infinity, while the coordinate z, for fixed r, 
increases monotonically at the interval (—00,00). The 
azimuthal angle tp ranges at the interval [Q, 2tt), as usual 
[22I ]. Also we assume that there exists at the spacetime 
an infinitesimally thin disk, located at the hypersurface 
z = 0, in such a way that the components of the metric 
tensor are symmetrical functions of z and that their first 
z-derivatives have a finite discontinuity at z = 0. 

Accordingly with the above considerations, 

9ab(r, -z) = g ab (r,z), (1) 
in such a way that, for z =/= 0, 

g a b,z{r, -z) = -g a b,z( r ' z )- ( 2 ) 
Thus then, the metric tensor is continuous at z = 0, 

[9ab] = 9ab\ z = Q+ - 9ab\ z = _ = 0, (3) 

whereas that the discontinuities in the derivatives of the 
metric tensor can be written as 

Jab = [&»&,*] = 2 9ab,z\ z=0+ , (4) 

where the reflection symmetry with respect to z = 
has been used. So, by using the distributional approach 
[H, IH, HI] , we can write the metric tensor as 

9ab = 9t b O(z)+g- b {l-0(z)}, (5) 

in such a way that the Ricci tensor can be written as 

Rab = R+ b e(z) + R- b {^ - 0{z)} + H ab S(z), (6) 

where 0{z) and 5{z) are, respectively, the Heaveside and 
Dirac distributions with support on z — 0. Here g^ b and 
iii are the metric tensors and the Ricci tensors of the 
z > and z < regions, respectively, whereas that 

Hab = \WI + YbSl - tfSffi - g zz j a b}, (7) 

where all the quantities are evaluated at z = + . 

Then, in agreement with ([5]), the energy- momentum 
tensor must be expressed as 

Tab = T+9(z) + T- b {l - 9(z)} + Q ab S(z), (8) 

where Tfr are the energy-momentum tensors for the z > 
and z < regions, respectively, and Q ab gives the part 
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of the energy-momentum tensor corresponding to the 
disk source. Accordingly, the Einstein equations, in ge- 
ometrized units such that c = 8irG — 1, are equivalent 
to the system 

Rab IdabR* = T±, (9) 
Hah — l^gabH = Q a b, (10) 

where H — g ab H a b and, again, all the quantities are eval- 
uated at z = + . Now then, when the thin disk is the 
only source of the gravitational field, so that Ti = 0, 
equation ([9]) reduces to the Einstein vacuum equations 

Rab = «, (11) 

for the z > and z < regions, respectively. Thus, in 
order to obtain solutions wtih a thin disk as source, we 
must solve the system ([TT]) by using in equation (fT0|) . 
as boundary conditions, the values of Q a b that describe 
properly the matter content of the disk. 

Now, in order to obtain explicit forms for the vacuum 
Einstein equations and the boundary conditions, we will 
take the metric tensor as given by the Weyl line element, 
written as [26] 

ds 2 = - e^dt 2 + e - 2 *[r 2 V + e 2A (dr 2 + dz 2 )}, (12) 

where $ and A are continuous functions of r and z. Fur- 
thermore, we will assume that $ and A are even functions 
of z, 

<$>(r,-z) = $(r,z), (13a) 

A(r,-z) = -A(r,z), (13b) 

in such a way that their first z-derivatives are odd func- 
tions of z, 

$ z (r,-z) = -$ z (r,z), (14a) 

A )Jt (r,-z) = -A z (r,z), (14b) 

which we shall require that do not vanish at z = 0. So, 
the vacuum Einstein equations (fTTj) are equivalent to the 
system 

(r*, r ) >r + (r$,,),« = 0, (15a) 
A, r = r($ 2 r -$ 2 2 ), (15b) 

A 2 = 2r$. r $ 2 , (15c) 

where (115a[) is the usual Laplace equation for an axially 
symmetric source in cylindrical coordinates, whereas that 
the intcgrability condition for the overdetermined system 
(I15bp - (|15cj) is granted when $ is a solution of (|15a|) , in 



such a way that A can be obtained by quadratures given 
a solution for 

On the other hand, equation (|10|) implies that the 
boundary conditions reduce to 

2e 2 (*- A > [A, - 2$ 2 ] = Ql (16a) 

2e 2 ^ A >A, 2 = Q*, (16b) 

where all the quantities are evaluated at z — + , and 
that Qab must have only two nonzero components. So, 
by using the orthonormal tetrad 

V a = e ~* <5 t a , (17a) 
W a = e* %/r, (17b) 

X a = e *-A ga^ (17c) 
ya = e *-A ( 1?d ) 

we can write the surface energy-momentum tensor Q ao 
in the canonical form 

Qab = eV a V b +pW a W b , (18) 

where e and p are, respectively, the energy density and 
the azimuthal pressure of the disk. In terms of these 
quantities, the boundary conditions can be written as 

2e 2 (*- A >[2$ z -A, 2 ] = e, (19) 

2e 2 ^- A >A, 2 = p. (20) 

Finally, by using (|15cp . we can cast these conditions as 
4e 2(*-A) [! _ r $ r j § z = £j (21) 

4e 2(*-A) r$ r$ 2 = P) (22) 

where, as before, all the quantities are evaluated at z — 
0+. 

As we can see from the above expressions, the more 
general energy-momentum tensor that is compatible with 
the line element (fT2j) and the boundary conditions (|10p . 
corresponds to a thin disklike source that only have a 
nonzero energy density and a nonzero azimuthal pres- 
sure. In agreement with this, instead of give specific 
prescriptions for the energy density and the azimuthal 
pressure, the Einstein equations will be solved only by 
requiring that these two quantities will be different from 
zero at the surface of a disk with an inner edge. Then, 
after a given solution be obtained, it can be used in order 
to obtain, from the boundary conditions, the correspond- 
ing expressions for the energy density and the azimuthal 
pressure. Therefore, the solution will correspond to the 
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more general static thin disk with an inner edge that can 
be obtained by exactly solving the Einstein equations. 

Accordingly, in order to obtain a solution representing 
a thin disk located in the hypersurface z = 0, with a 
circular central inner edge of radius a, we only need to 
imposse that 



*»(r,0+) 




(23) 



with f(r) any arbitrary function. Then, only after we 
find the more general solution, we will impose additional 
conditions in order to have a physically razonable behav- 
ior. So, in order to have an asymptotically flat spacetime, 
we will require that 



lim <3?(r, z) 

R — >oo 



0. 



lim A(r,z) = 0, 

R— >oo 



(24a) 



(24b) 



where R 2 



z 2 . Also, in order to have regularity at 



the symmetry axis, we will require that 
<&(0,z) < oo, 

A(0,z) < oo. 
We also will require that 

f(r) > 0, 

< r$ r < 1, 



(25a) 



(25b) 



(26a) 



(26b) 



in order that the energy density and the azimuthal pres- 
sure be positive everywhere. 

Finally, we need to check the finiteness of the total 
mass of the disks. So, in order to do this, first we take 
the mass density [i of the disks as defined by 



f = (Qat - \g ab Q)V a v\ 



(27) 



where Q = g ab Q a b and, as before, all the quantities are 
evaluated at z = + . Accordingly, by using (Til)]) and 
([15)). we have that the mass density reduces to 



fj, = e + p. 



(28) 



On the other hand, the total mass of the disks is given 
by 

/pItt poo 
/irfS = / / [j,e A - 2 *rdrd<p, (29) 
JO J a 

where dS = e A ~ 2 * rdrdip is the area element on the disk 
surface. So, by using (JSHJ, we can write the total mass 
as 



M = 2tt 



(e + p)e A - 2 *rdr, 



(30) 



III. SOLUTION OF THE EINSTEIN 
EQUATIONS 

In order to solve the Einstein vacuum equations, first 
we must solve the boundary value problem for $. How- 
ever, due to the nature of the boundary conditions fli"]). 
it is convenient to look for a different coordinate system 
that be naturally adapted to the geometry of the desired 
source. Accordingly, we introduce the oblate spheroidal 
coordinates as defined trough the relations 



= a 2 (l + a: 2 )(l-y 2 ), 



axy, 



(31a) 



(31b) 



where — oo < x < oo and < y < 1. So, when x = 
we have that z = and < r < a, whereas that when 
y = we have that z = and r > a. Furthermore, 
as x changes sign on crossing the surface y = 0, but 
do not changes in absolute value, this coordinate has a 
finite discontinuity when y = 0. Thus then, an even 
function of x is a continuous fuction everywhere but has 
a discontinuous normal derivative at y — 0. On the other 
hand, y is continuous everywhere. Therefore, the surface 
y = describes a thin disk with an inner edge of radius 
a, whereas that the surface x — describes the vacuum 
hole inside this edge. 

In the oblate spheroidal coordinates, the Weyl line el- 
ement (fT2l) can be rewritten as 



, (32) 



ds 2 = - e 2 *dt 2 + a 2 (l + x 2 )(l-y 2 )e- 2 *dy 2 

dx 2 dy 2 



+ a 2 (x 2 +y 2 )e 2(A -^ 



l + x 2 1-y 1 



in such a way that the Einstein vacuum equations reduce 
to 



[(l+.T 2 )$ x ], ;c + [(l- 2/ 2 )$,,],,=0, 



(33) 



the Laplace equation in oblate spheroidal coordinates, 
and the overdctcrmincd system 

A x = (l- y 2 )[x(l + x 2 )$%-x(l-y 2 )<Z> 2 v 

~2y(l + x 2 )<P iX $, y ]/(x 2 +y 2 ), (34) 

K y = (l + x 2 )[y(l + x 2 )<S> 2 x -y(l-y 2 )<P 2 y 

+ 2x(l~y 2 )<P^^ y ]/(x 2 +y 2 ), (35) 

whose integrability again is granted by equation (|33[) . 

On the other hand, by using (|31a[) and (|31bj) . is easy 
to see that 

$, x (0,y)/ay ; < r < a, 

(36) 

&, v (x, 0)/ax ; r > a. 



where the integration on if has been made. 



Accordingly, as the reflection symmetry of the solutions 
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implies that 



(37a) 



(37b) 



$ x (-^,y) = -&,x(x,y), 

the conditions (|2"3")) are equivalent to 

*,»(0,y) = 0, (38a) 



$ v (af,0) = F(a:); a: > 0, 



(38b) 



with -F'(x) any arbitrary function. The general solution 
of equation (f33|) with these boundary conditions is given 

by m 



= J^[iW*l»(l/) + B 2n Q2n{y)]p2n{x), (39) 
n=0 

where and i?2n are constants, i"2n(y) and Q2n(y) are 
the Legendre polynomials and the Legendre functions of 
the second kind, respectively, and pmix) = i~ 2n P2n{ix) . 
Therefore, all the solutions of the Einstein vacuum equa- 
tions for static spacetimes with any axially symmetric 
source as the considered here, are obtained by taking 
for the metric function <&(x,y) any particular choice of 
the above general solution, or expressions obtained from 
these solutions by means of linear operations. 

Now, in terms of the oblate spheroidal coordinates, 
condition (|24ap is written as 



lim 4>(x,y) = 0, 

x—*oo 

whereas condition (|25ap is written as 
$(x, 1) < 00. 



(40) 



(41) 



So, due to the behavior of the Legendre functions, it is 
clear that it is not possible to fulfill the physical con- 
ditions l|24p and (|25[) with any particular choice of the 
general solution (|39[) . However, by considering only the 
first term of the series, 



<f>o(x, y) = A + BoQo(y), 



(42) 



we obtain a solution that is regular for all y ^ 1. Then, 
if we take Aq = 0, this solution can be written as 



®o{x,y) = - In 



1 



y 



i-y 



(43) 



where a is an arbitrary constant, in such a way that a 
direct integration of (f3"4|) -([3"5 |) gives 



A o(x,y) = y hi 



i-y 2 



y 



(44) 



where the integration constant has been taken as zero. 
As we can see, this solution is not asymptotically flat 



neither regular at the symmetry axis. Nevertheless, by 
using (fTT)|) and we obtain for the energy density and 
the azimuthal pressure the expressions 



(4a/a)x 



2<>- 



P = o, 



(45) 



(46) 



in such a way that, if a > 0, the disk satisfy all the 
energy contions [28[. However, for any value of a 2 ^ |> 
the energy density increases without limit, may be at 
infinity or at the inner edge of the disk, whereas that 
for a 2 = | the energy density is everywhere constant, 
so that, in any case, the total mass of the disk will be 
infinite. 

On the other hand, although the previous solution has 
not a physically acceptable behavior, we can use it as 
the starting point to generate new and well behaved 
solutions. In order to do this, we consider the oblate 
spheroidal coordinates not only as functions of the cylin- 
drical coordinates (r, z), but also as parametrically de- 
pendents of the radius a, 



x = x(r,z;a), 



JJ = y(r,z;a). 



(47a) 



(47b) 



Accordingly, by considering also the metric function $ as 
dependent of a, 



<f> = $(r, z;a), 



(48) 



we can obtain a family of new solutions by applying the 
linear operation 



$„+i(r, z;a) 



d$ n (r,z;a) 
da 



(49) 



where n is an integer, n > 0. 

Thus then, by starting with the "seed solution" 
$>o(x,y), by means of the previous procedure it is gener- 
ated a family of new solutions that can be written in the 
simple form 

* n (r,z;a) =*„(,, y) = ^ffjffi^ , (50) 

for n > 1, where the F n (x,y) are polynomial functions, 
with highest degree An— 4, of which only we present below 
the first three, 

F, = 1, 



F 2 = x 4 + 3x 2 (l-y 2 )-y 2 , 

F 3 = 3x 6 (3-5y 2 ) + 5x 4 (6y 4 -lit/ 2 + 3) 
- x 2 y 2 (3y 4 - Sly 2 + 30) - y\y 2 - 3), 
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but all of them can be easily obtained by means of ((49 
So, is easy to see that 



and that 



lim <&„(£, y) = 0, 



$n(x,l) < 00, 



(51) 
(52) 



in fully agreement with conditions (|24aj) and (|25ap . 

Now, in order to obtain the correspondig metric func- 
tions A„(r, z; a), we make the integration 

A„(r, z;a) = A n (x,y) = J A, y (x,y)dy, (53) 

by taking A n (x, 1) = in order to grant regularity at the 
axis. So, by using (|5D|) in (|35[) . the obtained solutions 
can be written in the simple form 



A n (x,y) 



, -(2n-2)\(y 2 -l)A n (x,y) 
4 n a 2n (x 2 +y 2 ) 4n 



(54) 



for n > 1, where the A n (x,y) are polynomial functions, 
off highest degree 8n — 2, of which we present here only 
the first three, 

A, = x 4 (9y 2 -l) + 2x 2 y 2 (y 2 + 3) + y 4 (y 2 -l), 



A 2 = 2x L2 {9y 2 -I) -4x w (5ly 4 -Aly 2 + 2) 

+ x 8 (735y 6 - 1241/ + 419y 2 - 9) - x 6 y 2 (U2y 6 

- 1644?, 4 + 1604y 2 - 252) + x 4 y 4 {84y 6 - 384y 4 
+ 1266y 2 - 630) + Ax 2 y 6 (6y 6 + 6y 4 - 39y 2 + 63) 
+ 3yV + y 4 + y 2 -3), 

A 3 = 3x 16 (1225y 6 - 1275y 4 + 315y 2 - 9) 

- 24x 14 (980y 8 - 2095y 6 + 1205y 4 - 189y 2 + 3) 
+ 2x 12 (24255y 10 - 89475y 8 + 98472y 6 - 36316y 4 
+ 3473y 2 - 25) - 12x 10 y 2 (1835y 10 - 16665y 8 

+ 34716y 6 - 25292y 4 + 6001y 2 - 275) 

+ 6xV(90% 10 - 11946y 8 + 50563/ - 69397y 4 

+ 33365y 2 - 4125) + 8x 6 ?/ 6 (125y 10 + 926y 8 

- 9079y 6 + 24639y 4 - 22290y 2 + 5775) 

+ 6a;V(55y 10 + 29y 8 + 764y 6 - 4808y 4 + 8469y 2 

- 4125) + I2x 2 y 10 (5y 10 + 5y 8 + 80y 4 - 301y 2 

+ 275) + y 12 (5y 10 + 5y s + 5y 6 + y 4 + My 2 - 50), 

but all of them can be obtained as a result of compute 
the integral (|53[) . Accordingly, we have that 

lim A n (x,y) = 0, (55) 

x^oc 

and that 

A n (ai,l)=0, (56) 
in fully agreement with conditions (|24bp and (|25b[) . 



IV. BEHAVIOR OF THE SOLUTIONS 

As we can see from the expressions at previous sec- 
tion, by using (|49p was generated an infinite family of 
asymptotically flat solutions, all of them with a regular 
behavior at the symmetry axis. On the other hand, it 
seems to be that there is a singularity at the inner edge 
of the disks, when x = and y = 0. However, is easy 
to see that this apparent singularity is only a coordinate 
singularity. Indeed, as it is well known (26l . 29], the invari- 
ant characterization of the curvature must be in terms of 
14 scalars constructed from the Riemann tensor, R a bcd, 
and the metric tensor, g a b- So, in order to determine the 
nature of the aparent singularities, we must to compute 
all the curvature invariants for the family of solutions. 

Now, for any solution of the Einstein vacuum equa- 
tions, the curvature invariants are the 10 vanishing com- 
ponents of the Ricci tensor, R a b = 0, plus the four scalars 



v~ riabed / > 

K.I — K K a , 



K 



ii 



K 



in 



Kiv = 



bed > 

K fei-ft Kabcd, 

^ab ryklcd d 

6 klK Kabcd 



r ab r>kl rtmncd td 

£ kin. mn K Kabcd 



where g = det<? a ;, and e abcd is the Levi-Civita symbol. 
Furthermore, for any Weyl solution the last two invari- 
ants vanishes identically, so that we only need to compute 
K-j and fCn- 

By using the expressions (15"01) and ([54")) for & n (x, y) and 
A n (x, y), we can cast the curvature scalars as 



16a 2 e 4 ^- A ^N In (x,y) 



a 



(x 2 + y 2 ) 



2U2n 



(57a) 



K 



ii, 



48a 3 e e ^- A ^N IIn (x,y) 

„2U6n ' 



a 



8n+6( T 2 



(x 2 + y 2 y 



where Nj n {x, y) and Nn n (x, y) are polynomial functions, 
of highest degree 24n — 6 and 32rt — 9, respectively, that 
we do not write here explicitly due to their long size, but 
that all of them vanishes at the inner edge of the disks, 



Ar Jn (0,0)=AT //n (0,0)=0. 



(58) 



Morover, is easy to see that, in any neighborhood around 
(0,0), the difference between $> n {x, y) and A n (x, y) be- 
haves as 



- A r , 



a 2n (x 2 +y 2 ) 2n ' 



(59) 
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in such a way that 



lim 



= 4(<£>„-A„) 



(x l3 /)^(o,o) {x 2 + y 2 ) 12n 

p 6(<£>„-A„) 

lim 



0, 



0. 



(60a) 



(60b) 



(x, V )^(0,0) {x 2 + y2)16n 

and the limits exist, whatever be the path chosen to ap- 
proach the point (0,0). Accordingly, we have that 

(61a) 



lim K. In (x,y) = 0, 
lim JCun{x,y) = 0, 



(61b) 



and thus the curvature is regular at the inner edge of the 
disks. 

Now, in order to analyze the physical behavior of the 
sources, we will compute the energy density and the az- 
imuthal pressure for this family of disks. So, by using 
((3Ta|) and pTbj) . we can see that 



$ r (r,0) 



ax 



MM), 



r > a, 



and, by using ([50|) . is easy to prove that 
$ n , x (x,0)=0, n>l. 
Then, form ([22| . we can see that 
P« = 0. 



(62) 



(63) 



(64) 



That is, all the disks of the family have zero azimuthal 
pressure. 

On the other hand, by using equations (|2"Tj) . (|36| and 
(f50")l . the surface energy density of the disks can be writ- 
ten as 

a 2 (2n- 2)\B n {x)' 



4aE n (x) 

£n\X) — —-, — „ — 7T eXD 



(65) 



where x > and the E n (x) are positive definite polyno- 
mials of degree 2k, with fc = (n — l)/2 for odd n and 
k = n/2 for even n, of which we only will write below the 
first three, 

Ex{x) = 1, 
E 2 {x) = .x 2 + 3, 

E 3 (x) = 3(x 2 + 5), 

whereas that the B n (x) are positive definite polynomials 
of degree 4fc, with k = (n — l)/2 for odd n and fc = n/2 
for even n, the first three of them given by 

B 1 {x) = 1, 

B 2 {x) = 2x 4 + 8x 2 + 9, 
S 3 (£c) = 27a; 4 + 72x 2 + 50. 



From the above expressions we can see that, by taking 
a > 0, the energy density of the disks will be everywhere 
positive, 



,(*) > o. 



(66) 



So that, as the azimuthal pressure is zero, we have an 
infinite family of dust disks that are in fully agreement 
with all the energy conditions. Also is easy to see that, 
for any value of n, we have that 



,(0) 



lim e n (x) 



0. 



0. 



(67a) 



(67b) 



That is, the energy density of the disks is zero at their 
inner edge and vanishes at infinite. Furthermore, as the 
azimuthal pressure is zero, the mass density of the disks 
reduces to their energy density, 



Hn(x) = e n (x), 



(68) 



so that its behavior is the same as the energy density. 

Now, in order to show the behavior of the energy den- 
sities, we plot the dimensionless surface energy densities 
e n = ae„ as functions of the dimensionless radial coordi- 
nate f — r/a. So, in Figure [TJ we plot e„ as a function 
of f for the first three disks of the family, with n = 1, 2 
and 3, for different values of the parameter a n = a/a n . 
Then, for each value of n, we take a n = 0.5, 1, 1.5, 2, 
2.5, 3, 3.5 and 4. The first curve on left corresponds 
to a n = 0.5, whereas that the last curve on right cor- 
responds to a n = 4. As we can see, in all the cases 
the surface energy density is everywhere positive, having 
a maximun near the inner edge of the disks, and then 
rapidly decreasing as f increases. We can also see that, 
for a fixed value of n, as the value of a n increases, the 
value of the maximum diminishes and moves towards in- 
creasing values of r. The same behavior is observed for 
a fixed value of a n and increasing values of n. 

On the other hand, by using (|3T)|) . Q31ap and the 
total mass of the disks can be expressed as 



M n = 2ira 



f n {x)dx, 



where 



f n (x) = xe- A ^<°h n (x). 
So, fom (HU), d65j) and fTPl). is easy to see that 



,. fn+l(x) 
x ^°° Jn{X) 



tn+1 



(X) 



lim 

x^oo e n (x) 



= lim 



E n+ i(x) 



00 ax 2 E n (x) ' 

1 : n = 2k + 1. 



(69) 

(70) 

(71a) 
(71b) 

(71c) 







n = 2& + 2, 
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n= 1 




1 1.5 2 2.5 3 



n = 2 

3 | 1 1 



2.5 - 



2 - 



1.5 - 



0.5 - 








1 1.5 2 2.5 3 

n = 3 

1.2 | 1 1 1 1 



0.8 - 



C 




1.5 2 2.5 3 



FIG. 1: Surface energy density e n as a function of f for the 
first three disks of the family, with a„ = 0.5, 1, 1.5, 2, 2.5, 
3, 3.5 and 4. For each value of n, the first curve on left 
corresponds to a = 0.5, whereas that the last curve on right 
corresponds to a = 4. 

with k > 0. Therefore, by the limit comparison test for 
improper integrals, the convergence of M„+i is granted 
if M n is convergent and thus we only need to test the 
convergence of Mi . Indeed, a simple computation gives 

Mi = 2tt%/2^ T(l/4), (72) 

by granting so the convergence of all the mass integrals 
(155)1 . Accordingly, although the disks are of infinite ex- 
tension, all of them have finite mass. 



V. CONCLUDING REMARKS 

We presented an infinite family of asymptotically flat 
and everywhere regular exact solutions of the Einstein 
vacuum equations. These solutions describe an infinite 
family of thin dust disks with a central inner edge, whose 
energy densities are everywhere positive and well be- 
haved, in such a way that their energy-momentum ten- 
sor are in fully agreement with all the energy conditions. 
Moreover, although the disks are of infinite extension, all 
of them have finite mass. Now, as all the metric func- 
tions of the solutions are explicitly computed, these are 
the first fully integrated explicit exact solutions for such 
kind of thin disk sources. Furthermore, their relative 
simplicity when expressed in terms of oblate spheroidal 
coordinates, makes it very easy to study different dy- 
namical aspects, like the motion of particles inside and 
outside the disks and the stability of the orbits. 

Now, besides their importance as a new family of exact 
and explicit solutions of the Einstein vacuum equations, 
the main importance of this family of solutions is that 
they can be easily superposed with the Schwarzschild so- 
lution in order to describe binary systems composed by 
a thin disk surrounding a central black hole. Indeed, the 
superposition of the first member of this family with a 
Schwarzschild black hole already has been done, and was 
previously presented in [2l| , whereas that in a subsequent 
paper a detailed analysis of the corresponding superpo- 
sition for the full family will be presented. Accordingly, 
like was the family presented here, its superposition with 
the Schwarzschild black hole will be the first family of 
explictly integrated exact solutions for this superpostion 
of sources. 

Finally, it is worth to mention an interesting feature 
of this family of solutions. As in Newtonian theory the 
gravitational potential is given by the solution of the 
boundary value problem for the Laplace equation, we can 
consider the 3> n (r, z;a) as a family of Newtonian gravi- 
tational potentials of thin disklike sources with an inner 
edge, whose Newtonian mass densities are given by 



clearly diverging at the edge of the disks. Accordingly, we 
can conclude that there are not regular solutions within 
the Newtonian theory that properly describe the gravita- 
tional field of a thin disk with an inner edge, whereas that 
this kind of source can be properly described, by means 
of regular and asymptotically flat solutions, within the 
general relativistic garvitation theory 
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